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Abstract 

That /i ^ e,7 and r ^ /x,7 are sensitive probes of SUSY models with a see-saw mechanism is 
a well accepted fact. Here we propose a 'top-down' approach in a general SUSY 50(10) scheme. 
In this framework, we show that at least one of the neutrino Yukawa couplings is as large as the 
top Yukawa coupling. This leads to a strong enhancement of these leptonic flavour changing decay 
rates. We examine two 'extreme' cases, where the lepton mixing angles in the neutrino Yukawa 
couplings are either small (CKM-like) or large (PMNS-like). In these two cases, we quantify the 
sensitivity of leptonic radiative decays to the SUSY mass spectrum. In the PMNS case, we find 
that the ongoing experiments at the B-factories can completely probe the spectrum up to gaugino 
masses of 500 GeV (any tan (3). Even in the case of CKM-like mixings, large regions of the 
parameter space will be probed in the near future, making these two processes leading candidates 
for indirect SUSY searches. 
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I. INTRODUCTION 



It is known that Flavour Changing Neutral Current fFCNC) processes play an important 
role in the search for indirect signals of new physics [1]. On the other end the accumu- 
lating concordance between the Standard Model (SM) expectations and the vast range of 
experimental results in FCNC and CP violation point towards a low energy new physics 
which is flavour blind. For instance considering Supersymmetric (SUSY) extensions of the 
SM all this phenomenological evidence could be seen as a consequence of the fact that the 
mechanism that breaks SUSY and conveys SUSY breaking to the observable sector is com- 
pletely flavour blind. If this is indeed the case, it might be difficult to see any deviation at 
all from the SM predictions through FCNC j^. However, we already know that the SM is 
not enough in that it has to be supplemented by a mechanism to provide non-zero neutrino 
masses and mixings. An appealing example of this is the seesaw mechanism jsj where we 
introduce neutrino Yukawa couplings , and heavy right-handed neutrino Majorana masses 

Mr. 

In spite of the possible flavour-blindness of SUSY breaking, the Supersymmetrization of 
the Seesaw leads to new SUSY sources of Lepton Flavour Violation (LFV) ^ . The point is 
that in the running of the slepton mass matrices down to the electroweak scale, contributions 
to these matrices proportional to h'^h'^^ give rise to a mismatch in the diagonalisation of the 
lepton- slept on mass matrices and hence the appearance of flavour changing gaugino-lepton- 
slepton vertices. 

This effect was already noticed several years ago. What is new at this moment is the 
improvement in our knowledge of neutrino physics and prospects of better bounds from 
LFV physics. However, as shown in several recent studies , the low energy data on 

neutrino masses and mixings by themselves cannot predict the decay rates of Ij li'-f for 
a given set of SUSY breaking parameters. This ambiguity can be traced to the fact that 
the unknown parameters of h'^ and Mr cannot be completely fixed even after knowing all 
the three masses and three mixing angles of the neutrino sector. Thus any "bottom-up" 
approach to the study of these lepton fiavour violating processes suffers a large ambiguity 
which is encoded in an arbitrary orthogonal matrix R relating these low energy 
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neutrino parameters to the unknown Seesaw parameters and Mr^. Several interesting 
works Pi have studied this problem in detail for various neutrino spectra - hierarchical, 
inverse- hierarchical, degenerate - and under various assumptions on the mass eigenvalues of 
the right handed neutrino mass matrix. The bottom line of these analysis is that Ij —>■ li'-f and 
neutrino oscillation experiments provide complementary pieces of information to determine 
the seesaw parameter space jl^ . 

In this paper, we take an alternative point of view namely, we use a 'top-down' approach 
introducing a high energy framework which eliminates the ambiguity on the Seesaw param- 
eters. We consider a SUSY 5*0(10) model where all fermions in a generation are included in 
a single representation. The crucial find of our analysis is that in a generic SO (10) model 
at least one of the neutrino Yukawa couplings is of the same order as the large top Yukawa 
coupling. By generic we mean here that the result holds irrespective of the representations of 
the Higgs fields used to generate masses for the fermions. However, as expected the complete 
structure of h'^ remains unpredictable in such a general case. Essentially, the specification of 
the Grand Unified gauge group is not sufficient to uniquely determine the mixings appearing 
in the diagonalisation of h'^h"'' relevant for the analysis of Ij li,'y decays. 

Motivated by simple 5*0(10) models, we propose that the mixings diagonalising the com- 
bination h'^h'^'^ would have magnitudes within the two 'extreme' limits - namely, the CKM 
angles of quark mixing and the PMNS angles of the leptonic mixing. We show semi-complete 
50(10) models where these two 'extreme' cases can be realised. We then demonstrate that 
phenomenologically viable seesaw mechanisms can be realised in both these cases by a suit- 
able choice of the right handed neutrino Major ana mass matrix. We believe that these 
extreme cases would serve as "benchmark" scenarios for the seesaw induced lepton flavour 
violation within the context of SUSY 50(10). 

For both these "benchmark" scenarios, we have computed the Ij — >■ /j7 decay rates in a 
minimal supergravity (mSUGRA) or constrained MSSM (CMSSM) scenario, i.e, assuming 
completely universal soft SUSY breaking. We find that the present and future observational 
limits on — i> 67 can significantly constrain large regions in the parameter space even in 
the above mentioned small mixing (CKM) case. In the more optimistic case of large mixing 

^ It is possible that mechanisms other than the seesaw can be more predictive for Lepton Flavour Violation. 
See, for an example, Ref.Q. 
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(PMNS) the decay rates can become very high. Indeed, if the new proposals to explore LFV, 



mainly the /i — > 67 decay at PSI ll|, reach the planned sensitivity these LFV processes are 
going to be fully complementary with SUSY searches at LHC 

The rest of the paper is organised as follows. In the next section we show that in 5*0(10) 
models at least one of the neutrino Yukawa couplings is forced to be as large as the top 
Yukawa. Semi-complete SO{10) models are presented in section III for the minimal and 
maximal mixing scenarios. Results from numerical analysis and the constraints on SUSY 
parameter space are presented in section IV. Summarising remarks are presented in the last 
section. 



II. 50(10) AND NEUTRINO YUKAWA COUPLINGS 

As mentioned in the introduction, the top-down approach makes it possible to predict the 
neutrino Yukawa couplings removing the ambiguity associated with the seesaw parameters 
relevant for lepton flavour violation. In first place, we must choose a gauge group for the 
Grand Unified Theory (GUT). SO{10) is the natural choice as it is the minimal group that 
includes right handed neutrinos besides the rest of the SM fermions in a single representation. 
One would then have to consider a particular choice of the Higgs fields representations to 
generate fermion masses, which again introduces some sort of an 'ambiguity' at the high 
scale. To keep the discussion as general as possible, we do not resort to any particular 
model of fermion masses within SO{10), but try to see how well one can predict h'^ in a 
generic scenario. 

In the SO{10) gauge theory, all the known fermions and the right handed neutrinos are 
unified in a single representation of the gauge group, the 16. To analyse the Yukawa matrices 
in this framework, we need to specify the superpotential. In principle, the superpotential 
can receive contributions both from renormalisable and non-renormalisable terms The 
product of two 16 matter representations can only couple to 10, 120 or 126 representations 
which can be formed either by a single Higgs field representation or a non-renormalisable 
product of representations of several Higgs fields. In either case, the Yukawa matrices 
resulting from the couplings to 10 and 126 are complex symmetric whereas they are anti- 
symmetric when the couplings are to the 120. 

Therefore, the most general 5*0(10) superpotential relevant for fermion masses can be 
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written as 

Wsom = h]!ll6i 16j 10 + h]fH6i 16, 126 + h]fl6i 16, 120, (1) 

where i,j refer to the generation indices. In terms of the SM fields, the Yukawa couphngs 
relevant for fermion masses are given by jlS^ : 

16 16 10 D 5 (mm' + z/z/') + 5 {dd'' + ee"), 
16 16 126 D 1 z/'z/' + 15 z/z/ + 5 (m' - 3 uu^) + 45 {dd" - 3 ee"), 
16 16 120 D 5 z^z/" + 45 mm' + 5 {dd" + ee') + 45 {dd" - 3 ee'), (2) 

where we have specified the corresponding SU{5) Higgs representations for each of the 
couplings and all the fermions are left handed fields. From above, it is clear that if only 
the 10 and 126 Higgs representations are present in the theory, the Yukawa matrices of the 
down quarks and charged leptons as well as the up quarks and neutrinos are deeply related. 
In fact, in a model where only the 10-plets are present we would have exact quark-lepton 
Yukawa unification, not only among charged leptons and down quarks (6-r unification) 
but also for up quarks and (Dirac) neutrinos. Similarly a dominant contribution from the 
126 representation would predict as well quark-lepton unification, although introducing the 
Georgi-Jarlskog factors of 3 If only one of these representations or any combination 
of them contributes, the Yukawa matrices would be exactly symmetric. These properties 
are broken by the introduction of the 120 representation. The Yukawa couplings of this 
representation are anti-symmetric and break the quark-lepton unification because they can 
contribute independently to the quark and lepton Yukawa matrices. In general, both the 
symmetric and anti-symmetric contributions can be present leading to Yukawa matrices of 
generic nature. 

The resulting mass matrices can be written as 

= Ml + Mf26 + Mflo, (3) 

Mln = Ml-3Ml, + Ml,, (4) 

= M^o + + Mf,, + Mg, (5) 

M' = Mfo - 3Mg + m!,, - 3Mfl„ (6) 

ML = Mil,, (7) 

M^n = Ml,,. (8) 
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We come now to discuss the main result of this section: At least one of the Yukawa 
couplings in h" = ^lr large as the top Yukawa coupling. This resuh holds 

true in general independently from the choice of the Higgses responsible for the masses in 
Eqs. (jni El) provided that no accidental fine tuned cancellations of the different contributions 
in Eq. (0)) are present. If contributions from the lO's solely dominate, h'^ and /i" would 
be equal. If this occurs for the 126's, then h'^ = —3 h^. In case both of them have 
dominant entries, barring a rather precisely fine tuned cancellation between M^q and Mf2Q 
in Eq. Q, we expect at least one large entry to be present in h'^. A dominant antisymmetric 
contribution to top quark mass due to the 120 Higgs is phenomenologically excluded since 
it would lead to at least a pair of heavy degenerate up quarks. However, the 120 can still 
provide a non negligible contribution to the up quark masses and in particular, to that of 
the top quark, if there is at the same time a large symmetric contribution. In this case 
the up quark and neutrino Yukawa matrices have both large symmetric and antisymmetric 
contributions. The above stated result holds also in this general situation (for a complete 
proof, see the Appendix). 

Apart from sharing the property that at least one eigenvalue of both M" and M]^^ 
has to be large, for the rest it is clear from Q and ^ that these two matrices are not 
aligned in general, and hence we may expect different mixing angles appearing from their 
diagonalisation. This freedom is removed if one sticks to particularly simple choices of the 
Higgses responsible for up quark and neutrino masses. For instance, as long as one sticks to 
representations of the Higgs fields which preserve the quark-lepton Yukawa unification, like 
the 10-plets (due to the underlying Pati-Salam symmetry), the mixing angles appearing in 
the diagonalisation of the up-quark mass matrix also appear in h'^. In this case h'^ can be 
completely predicted at the high scale. However, in a general scenario, this need not hold 
true. Large contributions from 120 can break this alignment between up and the neutrino 
Yukawa matrices as is evident from Eqs. (jHl-El) and hence h'^ would be no longer predictable 
in this general situation. 

Keeping with our philosophy to be as general as possible, we find two cases which would 
serve as 'benchmark' scenarios for seesaw induced lepton fiavour violation in SUSY 5*0(10). 
The first one corresponds to a case where the mixing present in h'^ is small and CKM-like. 
This is typical of the models where fermions attain their masses through 10-plets. We will 
call this case, 'the minimal case'. As a second case, we consider scenarios where the mixing 
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in h'^ is no longer small, but large like the observed PMNS mixing. We will call this case 
the 'the maximal case'. Within SO{10) this is possible in models with asymmetric Yukawa 
matrices. These two 'benchmark' cases, we believe, would span the range of lepton flavour 
violation generated by h'^h'^'^ at the weak scale in a fairly general way. We now proceed to 
study these constraints for the two cases. 



III. 50(10) AND LEPTON FLAVOUR VIOLATION 

As mentioned in the Introduction, the off-diagonal entries in the slepton mass matrices 
are generated by the product h'^h'^ through renormalisation group evolution^. Having 
known the matrix h'^ or even the product h'^h'^ ^ at the high scale, the off-diagonal entry 
[Tn'jjij can be calculated at the weak scale. These entries would then contribute to the lepton 
flavour violatmg decays I, ^ l.,r A naive estimate of the branehiug ratios of these deeays 
is given by (in standard notation) [15j: 

BR(/,-/.7)-4^^^#^tan^A (9) 
where msusY represents the typical soft Supersymmetric breaking mass. The present ex- 



perimental limits on the branching ratios of these decays are given as 



16 



III 



Si?(/i e, 7) < 1.2 X 10"^\ (10) 
BR{t /i,7) < 5. X 10"^ (11) 

In future these bounds are expected to improve at least by a few orders of magnitude. In 



particular, in t 
to improve to 



le proposed experiment at PSI, the limits on the BR{fi e, 7) are expected 
HI 

BR{fi ^ e,-f) <10-^\ (12) 



These limits would now constrain the parameters governing the sparticle spectrum, namely 
mo, Ml/2, ^0, sg(/i) and tan/5 in the mSUGRA with electroweak radiative breaking. It should 
be noted that the above mentioned RG effects would require the scale of Supersymmetry 
breaking to be higher than the scale of right handed neutrinos. Thus in models of Gauge 



^ We have implicitly assumed mass matrices to be real for the following discussion and thus neglect any 
effects due to the presence of phases in the mass matrices. 
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Mediated Supersymmetry Breaking (GMSB) these effects would be absent. We now analyse 
these branching ratios in detail for the two cases mentioned at the end of the previous 
section, namely, the small mixing case with CKM angles and the large mixing case with 
PMNS angles in h". 



A. The minimal Case: CKM mixings in 

The minimal Higgs spectrum to obtain phenomenologically viable mass matrices includes 
two 10-plets, one coupling to the up-sector and the other to the down-sector. In this way 
it is possible to obtain the required CKM mixing in the quark sector. The 5*0(10) 
superpotential is now given by 

Wsom = \ KflQ^ 16, 10„ + ^ }4fl&, 16, 10, + ^ /ig 16, 16, 126. (13) 

We further assume the 126 dimensional Higgs field gives Majorana mass only to the right 
handed neutrinos. An additional feature of the above mass matrices is that all of them 
are symmetric. Without loss of generality we can rotate the 16-plet into a basis where the 
charged leptons and the down-type quarks are diagonal. In terms of the SM fields, we can 
rewrite the above as 

W = Kt^ Q, u'j Hu + hi Q, d1 Hd + hi el + h\, U v] Hu + ^M^,, vl v]. (14) 

Immediately we see that the following mass relations hold between the quark and leptonic 
mass matrices at the GUT scale^: 

/i" = /i^ ; h'^ = h\ (15) 
In the above basis, the symmetric h^ is diagonalised by: 

^CKM /i" ycKM = Kiag- 

Hence from (jl5j) : 



VcKM /i" VcKM - Kiaq- (16) 



h'' - VcKM Kiag VcKM- (17) 



^ Clearly this relation cannot hold for the first two generations of down quarks and charged leptons. As 
usual, small corrections due to non-renormalisable operators or suppressed renormalisable operators |l4| 
can be invoked. 
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According to Eq. (0), BR(/i — > 67) depends on: 

[h''h'']2^^h'',Vt,Vts + Oihl). (18) 

In this expression, the CKM angles are small but the presence of the large top Yukawa 
coupling compensates for such suppression. The large couplings in h'^ ~ 0{ht) induce 
significant off-diagonal entries in m| through the RG evolution between Mqut and the scale 
of the right-handed Majorana neutrinos ^, Mr.. The induced off-diagonal entry relevant for 
yU — > e, 7 is of the order, 

K]2i ^ -^(3^0 + ^0) hi Vu yts\og{^^) + 0{hl). (19) 

The required right handed neutrino Majorana mass matrix consistent with both the observed 
low energy neutrino masses and mixings as well as with CKM like mixings in h'^ is determined 
easily from the seesaw formula defined at the scale of right handed neutrinos as 

= -K" ^ vl (20) 

= -h^ h'^ vl (21) 

where we have used the symmetric nature of the h"^ in the second equation. Inverting 
Eq. (j2ni), one gets: 

Mr = -h" m-' vl 

= VcKM Kiag VcKM f^u^VcKM Kiag ^CKM^ (22) 

where we have used Eq. fll7|l for h'^ . Furthermore, can be written as = 

UpMNS diag[ni~^] UpMNSi whose entries are determined at the low scale from neutrino 
oscillation experiments. The structure of Mr can now be derived^ for a given set of neutrino 
masses and mixing angles. Neglecting the small CKM mixing in h'^ we have 

^ K[^v\i Khc[m~^]i2 Kht[m^^]i;i'^ 



(23) 



^ Typically one has different mass scales associated with different right handed neutrino masses. 

^ The neutrino masses and mixings here are defined at Mqut- Radiative corrections can significantly 
modify the neutrino spectrum at the weak scale This is more true for the degenerate spectrum of 
neutrino masses ji^l and for some specific forms of For our present discussion, with hierarchical 

neutrino masses and up-quark like neutrino Yukawa matrices, we expect these effects not to play a very 
significant role. 
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It is clear from above that the hierarchy in the Mr mass matrix goes as the square of the 
hierarchy in the up-type quark mass matrix. Furthermore, for a hierarchical neutrino mass 



UpMNS 



(24) 



spectrum we have rriy^ ^ yAm^^^, m^^ ~ v^"^© ^^"^ ^j^i ^ y^"^© for a nearly 
bi-maximal Upmns '■ 

1/V2 \ 

-1/2 1/2 
V 1/2 -1/2 I/V2/ 
it straight-forward to check that all the right handed neutrino mass eigenvalues are controlled 
by the smallest left-handed neutrino mass. 

M^3 - ; ^ ; Mp,^^ . (25) 

4 rriy^ 4 m,,^ 2 m,,^ 

This implies that we can not choose an arbitrarily small neutrino mass if we want the right- 
handed neutrino masses to be below Mqut- In our numerical examples in section IV, we 
choose mjyg = 0.05 eY^m^^ = 0.0055 eV, m,^j = 0.001 eV. 

At this point we have both h"^ and Mp determined and we can now use the experimental 
bounds on BR(/i e, 7) to constrain the SUSY parameter space. The only unknowns are 
the SUSY breaking soft terms mo, M1/2, Aq and sg(/i), tan/5. For instance, from Eqs. ()9I19|) 
we have 

"""-Stt^G^ '^^^M^' B ^^^^ 
where B represents the experimental limit on the branching ratio BR{^ ^,7), Aq is 
assumed to vanish and mo is identified with msusY- Taking the futuristic limit of B < 10~^^, 
we see that mo can be probed up to 1 TeV for large tan /? ~ 40, bordering the limit that 
will be probed at LHC [22]. We quantify these results with a numerical analysis in the next 
section. 



B. The maximal case: PMNS mixing angles in 

The minimal 5*0(10) model presented in the previous sub-section would inevitably lead to 
small mixing in h'-'. In fact, with two Higgs fields in symmetric representations, giving masses 
to the up-sector and the down-sector separately, it would be difficult to avoid the small CKM 
like mixing in h'^ . To generate mixing angles larger than CKM angles, asymmetric mass 
matrices have to be considered. In general, it is sufficient to introduce asymmetric textures 
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either in the up-sector or in the down-sector. In the present case, we assume that the down- 
sector couples to a combination of Higgs representations (symmetric and anti-symmetric) ^ 
leading to an asymmetric mass matrix in the basis where the up-sector is diagonal. We 
have : 

Wsom = \ Kt 16. 16,10" + ^ h!lf 16, 16,$ + i /i« 16, 16,126 , (27) 



where the 126, as before, generates only the right handed neutrino mass matrix. To study 
the consequences of these assumptions, we see that at the level of S't/(5), we have 



Wsu{^) = \ hi 10. lOi 5„ + K 5i U 5u + h% 10, 5,- 5^ + ^ 1,1^, (28) 



where we have decomposed the 16 into 10 + 5 + 1 and 5.. and are components of 10„ and 
$ respectively. To have large mixing ~ Upmns in h'^ we see that the asymmetric matrix 
h'^ should now be able to generate both the CKM mixing as well as PMNS mixing. This is 
possible if 

^CKM UpMNS — hdiag- (29) 

This would mean that the 10 which contains the left handed down-quarks would be rotated 
by the CKM matrix whereas the 5 which contains the left handed charged leptons would be 
rotated by the Upmns matrix to go into their respective mass bases [23]. Thus we have, in 
analogy with the previous sub-section, the following relations hold true in the basis where 
charged leptons and down quarks are diagonal: 

h"" = VcKM Kiag ^CKM ) (30) 

h'^ = UpMNs hl,^. (31) 
Using the seesaw formula of Eq. and Eq. we have 

M, = D,ag{^. (32) 

mi,^ rriy^ m^.j^ 

This would mean that this setup would require Mr to be diagonal at the 5*0(10) level in 
the basis of diagonal K^'^ , Eq. ()27|) . We now turn our attention to lepton flavour violation 
in the scenario. The branching ratio, BR(yU e, 7) would now be dependent on: 

[h-W ^]2i = h] U,s + hi U,2 f/e2 + Oihl), (33) 



^ The couplings of $ in the superpotential can be either renormaHsable or non-renormaUsable. See |23| for 
a non-renormalisable example. 

11 



where Ufi are elements of the Upmns matrix. It is immediately clear from the above that 
in contrast to the CKM case here the dominant contribution to the off-diagonal entries 
depends on the unknown magnitude of the element Ues- If Ues is very close to its present 



limit ~ 0.2 2^, the first term on the RHS of the Eq. (j33|) would dominate. Moreover, this 
would lead to large contributions to the off-diagonal entries in the slepton masses with 
of C(l). We have 

[mlhi ^ -^(3m^ + Al) hi ^e3log(^) + 0{hl). (34) 

The above contribution is large by a factor {Uf^-^Ue^) / {YtdVts) ~ 140 compared to the CKM 
case. From Eq. Q we see that it would mean about a factor 10^ times larger than the CKM 
case in BR(/i 6,7). In case t/es is very small, i.e, either zero or ~ h^/hf Ue2 ~ 4 x 10~^, 
the second term oc h"^ in Eq. (j33|) would dominate. However the off-diagonal contribution 
in slepton masses, now being proportional to charm Yukawa could be much smaller, in fact, 
even smaller than the CKM contribution by a factor 

If f/e3 is close to it's present limit, the current bound on BR(/i — 6,7) would already be 
sufficient to produce stringent limits on the SUSY mass spectrum. Indeed from: 

"^"^S^l^^^il^" ^^tan/5 (Ao = 0), (36) 
B < 10^^^ probes tuq at the TeV level even for small tan/5. We make these statements 
more concrete in the next section with our results from numerical analysis. 



IV. 50(10) AND SUSY MASS SPECTRUM 

As mentioned in the introduction we have chosen to work within a mSUGRA framework 
with flavour blind universal soft breaking terms at the GUT scale. It is known 12] that if 
the soft breaking terms are hard at scales above Mqut and if the perturbative RG approach 
can be safely used in such energy interval above Mquti then the large value of ht can 
lead to relevant radiative contributions to the mass of the stau, spoiling the slepton mass 
universality at Mqut- Here we are interested in the effects produced by the large h'^ coupling 
in the running of the m| down to the Mr scales which are below Mqut- For simplicity, 
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we assume the soft masses to be universal at Mqut, hence neglecting the possible above 
mentioned effects which would further enhance the stringent bounds that we obtained here. 

In our numerical analysis, we have considered MSSMRN (MSSM + Right handed neu- 
trinos) RGE from the scale Mqut down to the scale of the different RH neutrinos, which 
we integrate out in several steps until the scale of the lightest RH neutrino, M^^. /^From 
Mr^ to Mz the standard MSSM RGE are used ^|. As a result of these renormalisation 
group effects, the slepton mass matrices which were diagonal at Mqut to start with are 
now non-diagonal at the weak scale. At this latter scale, we have numerically diagonalised 
these mass matrices and the corresponding eigenvalues and mixing matrices are found. We 
used the complete calculations of Hisano et. al |3| to compute the branching ratios of both 
/i — >• e, 7 and r yU, 7. 

We produce scatter plots of BR(/i — »• e, 7) and BR(r — »• /x, 7) vs. M1/2 for small and large 
values of tan (3. We restrict the allowed SUSY parameter space by imposing i) experimental 
constraints from direct sparticle searches, in particular requiring m^- and m^± to be above 
Mz and ruu > 33 GeV; ii) the LSP to be neutral and iii) 2 x 10"^ < BR(6 ^5,7) < 4 x 10"^. 
In each scatter plot, for a given tan f3 and sg{fi), the rest of the parameters are allowed to 
vary within the ranges: 

a) 90 GeV < mo < 900 GeV ; 

b) 90 GeV < Ml/2 < 700 GeV, 

and c) Aq is parameterised as Aq = Ai ttiq with Ai lying between (-3, 3). For the neutrino 
mass eigenvalues, we have used the values specified at the end of section IIIA. We have 
chosen sin^ ^ 0.71; sin^ 26'atni ~ 1 ^^id Ue3 ~ 0.15. 

In Fig. la) and lb) we show the scatter plots for BR(/i — e, 7 ) for the CKM case and 
tan (3 = 2 and tan /5 = 40 respectively. Similar figures for BR(r yU, 7 ) and the same 
values of tan/? are presented in Figs. Ic) and Id). All these plots are calculated with /i > 
but the results do not change significantly with negative /x. These plots reflect an interesting 
correlation between the branching ratios and the GUT value of the universal gaugino mass. 
This is due to the fact that the gaugino mass fixes the chargino and neutralino masses at 
M\Y and, to a small extent it also influences the slepton masses through RGE. However, 
for a fixed M1/2 the different values of mo and Aq can change the value of the BR within a 
range of 3 orders of magnitude. Nevertheless, this fact has still important consequences. For 
instance, for tan/3 = 40 reaching a sensitivity of 10"^'^ for BR(/i 67) would allow us to 
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probe completely the SUSY spectrum up to M1/2 = 300 GeV (notice that this corresponds 
to gluino and squark masses of order 750 GeV) and would still probe a large regions in 
parameter space up to M1/2 = 700 GeV. In the case of smaller values of tan/5 the BR 
scales as (tan/5)^ and therefore for tan/5 = 2 only a small part of the parameter space with 
Ml/2 < 300 GeV can be probed with this sensitivity. Similarly in the r — /X7 decay a 
sensitivity of 6 x 10~^ which could be reached in the B factories in the near future [itI I25I1. 
would allow to probe a sizeable piece of the parameter space for large tan (3. 

In the PMNS scenario Figs. 2a) and 2b) show the plots for BR(/i e, 7) for tan (3 = 
2, 40, whereas Figs. 2c and 2d are for BR(r /i,7). As we said in the previous section 
(see Eq. (jSlI)) in the PMNS case, the results concerning BR(/i — > 6,7) strongly depend on 
the unknown value of Ues- In the plots 2a and 2b, the value of Ues chosen is very close to 



the present experimental upper limit j2^. As long as Ues ~ 4 x 10"^, the plots scale as U^^, 
while for Ues ~ 4 x 10~^ the term proportional to in Eq. ()34j) starts dominating and 
then, the result is insensitive to the choice of t/es- Here we see that with the present limit 
on BR(yU 6,7), all the parameter space would be completely excluded up to M1/2 = 700 
GeV for Ues = 0.15 for any value of tan/?. For instance, a value of Ues = 0.01 would reduce 
the BR by a factor of 225 and still most of the parameter space for tan/5 = 40 would be 
completely excluded. For tan/5 = 2 this would probe approximately half the parameter 
space up to M1/2 = 300 GeV. Contrary to expectations the present bound of the t fij 
starts exploring the SUSY parameter already for low tan/5 = 2. For larger tan /5 = 40 even 
the present bound rules out significant regions of the parameter space. The main advantage 
of this decay mode is that it does not depend on the value of f/es and therefore provides an 
important constraint on the parameter space of the model for any value of t/es- 

It is important to emphasise that these radiative leptonic decays are a much more powerful 
probe on the SUSY parameter space than the very similar 6 — 57 decay. This can be traced 
to the fact that slepton and chargino/neutralino masses are a factor -\/6 smaller that the 
corresponding gluino and squark masses for the same GUT initial parameters due to RGE 
effects. Therefore, these decays will become in the near future the most stringent constraints 
on the SUSY parameter space and offer an excellent opportunity for SUSY searches. 
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V. DISCUSSION OF THE RESULTS AND FINAL REMARKS 



As we have mentioned, 5*0(10) is an interesting and predictive scheme for BR(/j hjl) 
in so that one of the neutrino Yukawa couphngs has to be of the 0{ht). In our plots of Figs 
1 and 2 we exhibited resuhs of the two 'extreme' cases of "small " (CKM-like) and "large " 
(PMNS -like) lepton mixing angles. We see that the constraints coming from BR(/i e, 7) 
and BR(r Ai, 7) are very significant. In particular for the PMNS case, even the present 
bounds on these two branching ratios (Eqs. (HnHHI)) are able to exclude large regions of the 
SUSY 5*0(10) parameter space which are still allowed by all the present SUSY accelerator 
tests. 

Needless to say, improving the experimental sensitivity of those two processes, would 
place lepton radiative decays in the forefront of the indirect searches of SUSY signals. In- 
deed, reaching O(10~^^) for BR(/i 6,7) would probe larger regions of the SUSY 5*0(10) 
parameter space than that probed by 6 — s, 7 even in the less optimistic case of the small 
CKM lepton angles. Obviously, this statement becomes stronger when one moves to the 
PMNS case unless Ues is really very small, say f/es ~ 10~^. If this latter circumstance oc- 
curs, then T becomes more powerful SUSY probe than /i — e, 7 and again, reaching 
sensitivity for BR(r /i,7) of C(10^^) we could do better than 6 — >• 5,7 in constraining 
the SUSY 50(10) parameter space. 

Our analysis, referring in all generality to the framework of a relevant class of unified 
SUSY models, once again emphasises the importance and need for a strenuous effort in 
pursuing the challenging experimental road of lepton radiative decays. 
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Appendix I 

In this appendix we discuss in detail the proof that the large hierarchy in M" necessarily 
requires that the symmetric Yukawa couplings include a large entry of 0{mt). Let us 
decompose a generic mass matrix in symmetric and anti-symmetric parts as follows: 

M = Ms + Ma, (37) 
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where 



Mo 



( an ai2 ^ 

ai2 (122 ^23 

y ai3 a23 033 / 



and 



-bi2 623 

V -bl3 -&23 J 

The eigenvalues of these matrix can be found solving the cubic equation: 







bi2 bi3 ^ 



X-^ + aX^ + bX + c^O 



with the coefficients, 



a = Tr [M"M" ^] 

= aj^ + + + 2 {a\^ + 0^3 + 0^3 + b\^ + 6^3 + b\^, 



(38) 



(39) 



(40) 



(41) 



(42) 



where Aljj[M"M" ^] is the minor obtained by omitting the % row and j column from the 
determinant. These are given as 

Mxx = (^22 + (^12 - 012)^ + (023 + h^-if) {{bi3 - anf + 0^3 + (623 - ^23)^) 

- (022 (023 - ^23) + (023 + ^'23) 033 + (612 - 012) (&13 - 013))^ (43) 

A^22 = (an + (&12 + 012)^ + (ai3 + ^13)^) ((&13 - a^f + 033 + (^23 - ^23)^) 

- (an (ai3 - &13) + (^13 + ^'13) 033 + (612 + 012) (023 - ^23))^ (44) 



M33 = (an + {bi2 + auf + (ai3 + 613)' ) ((&12 - 012)' + 4^ + (623 + 023)') 

- (an (ai2 - &12) + (012 + &12) ^12 + (^13 + &i3)(a23 + ^'23))^ (45) 

The parameter c is given as follows: 

c = Det [M"M" ^] 

~ (&23 0^11 ~ 2 bi3 623 0^12 + 2 bi2 &23 ^^13 + ^13 ^^22 ~ ^^3 022 ~ 2 612 613 023 

+ 2 ai2 ai3 a23 - Oii ^23 + ^12 «33 - 012 «33 + On ^22 033) (46) 
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To generate a hierarchical spectrum, as required by the up-quark sector, a, b and c have to 
satisfy the following conditions: 



2 I 2 I 2 

a = m„ + + 

r, 2 2,2 2,2 2 

2 2 2 

c = m^m^m^ 



(47) 
(48) 
(49) 



We want to prove that at least one of the elements in the symmetric matric must be 
of order m^. To do this we show that the case where rrit has its origins solely from the 
anti-symmetric part i.e, the elements bij is not consistent with the observed spectrum. In 
such a case at least one of the bij is as large as However from Eqs. ()41ll46p we see that 
such an assumption would be in conflict with the condition in Eq. ()48|) as it leads to 6 ~ mf 
through one of the A4ij. For ex: if fo^g ~ mf ; Ain ~ mf. This is in fact true for any of 
the bij. Therefore a large bij would lead to a degenerate spectrum rather than a hierarchical 
spectrum required. Thus rrit cannot have its origins solely from the anti-symmetric part. 
The only allowed possibilities would then be, either rrit has its origins solely in the symmetric 
part or through elements of both Ms and M4, when some elements in both matrices are of 
comparable magnitude ~ 0{mt). In either case, it is clear that Ms would contain at-least 
one element as large as mt- 
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FIG. 1: The scatter plots of branching ratios of rare leptonic radiative decays vs. M1/2 are shown 
for the (minimal) CKM case for two specific values of tan (3. Results do not alter significantly with 
the change of sign(;^). 
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FIG. 2: The scatter plots of branching ratios of rare leptonic radiative decays vs. M1/2 are shown 
for the (maximal) PMNS case for two specific values of tan f3. Results do not alter significantly 
with the change of sign(/i). 
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